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Abstract 

We have proposed a new method for the Coulomb wave function correction 
which includes the momentum resolution for charged hadron pairs and applied 
it to the precise data on tt +/ k~ correlations obtained in p + Ta reaction at 
70 GeV/c. It is found that interaction regions of this reaction (assuming 
Gaussian source function) are 9.8 ±5.8 fm and 7.7 ±4.8 fm for the thicknesses 
of the target 8 /im and 1.4 /im, respectively. We have also analyzed the data 
by the formula including strong interaction effect. The physical picture of the 
source size obtained in this way is discussed. 
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I. INTRODUCTION 



We have obtained recently new formulae for the Coulomb wave function corrections for 
charged hadron pairs |l],j2|. Among other things we have applied them to data on 7r + 7r~ 
correlation obtained in p + Ta reaction at 70 GeV/c (and corrected by the usual Gamow 
factor only). However, as was pointed out to us by one of the author ||, our approach does 
not account for the finite momentum resolution published in |J. In fact, our formulae cannot 
be applied directly to experimental data in which such momentum resolution is accounted 
for. Therefore in the present work we would like to extend our method for the Coulomb 
wave correction provided in [T],Q] in such way as to include also the momentum resolution 
case and we would like to re-analyse data of as well as the new data of || obtained with 
two kinds of thickness of the Ta target: 8 /xm and 1.4 /xm. 

In the next Section we provide, for the sake of completeness, the main points of the 
analysis method used in Ref. || whereas in the Section 3 we outline our new method for the 
Coulomb wave correction including this time also the momentum resolution. In Section 4 
we consider also the effect of the strong final state interactions. The final part contains our 
concluding remarks. 



II. DATA RECONSTRUCTION METHOD (USING THE GAMOW FACTOR 

WITH MOMENTUM RESOLUTION) 

It has been stressed in || that relative momenta of tt + tt~ pairs observed by them have 
some finite resolutions. The averaged correlation function, defined as 



*(*+,-) = / - - = G„l + (1 - !)•(,), (1) 




depends therefore on this momentum resolution, where G c denotes a correlation function 
of 7r + 7r~ pairs from decays of short-lived resonances (SLR), and $(?) denotes a correlation 
function between pions from decays of the long-lived resonances (LLR) or LLR's and SLR's. 
The parameter b means the portion of the correlation function governed by the Coulomb 
wave function. For the sake of simplicity $(<?) = 1. Fig. |I](a) shows physical picture of Eq. 
(H|). To account for it the following random number method has been proposed in Ref. || 
in order to obtain the corresponding averaged quantities in analysis of the correlation data: 

(i) The relative momentum of the measured pair, q = pi — p 2 , is decomposed into its 
longitudinal and transverse components: q L and q T , respectively, by making use of the 
uniform random number u G (0, 1) (it is worthwhile to notice at this point that the 
transverse components q T in data of are smaller than 10 MeV/c). The following 
scheme has been employed: 

J 10^ for g>10MeV/c, 
qT ~\qV^ for q < 10MeV/c , 1 a) 

q, = Jq 2 - ql ■ (2b) 
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[ii] The Gaussian random numbers for g L and q T are generated in the following way: 



Ql (random) X + g L , (3) 

<?T(random) = &tX 4~ Q T . (4) 

where X stands for the standard Gaussian random number || whereas a L and a T are 
longitudinal and transverse setup resolutions for the corresponding components, which 
are equal to (values used in ||): a L = 1.3 MeV/c; <r T = 0.6 MeV/c (for target of the 
thickness 8 /im ) and cr T = 0.4 MeV/c (for the 1.4 pm target). 



(Hi) Using the randomized number g ran dom = y ^(random) + ^(random) one calculates tne 
corresponding randomized Gamow factor correction: 

(-it \ 27r^ ranc jom ir\ 

exp(-27n? random ) - 1 



where Random = ma/q random . 

Calculating now the average value of G(— ^random) in 10 ~ 100 k events one can estimate the 
Gamow factor with this finite momentum resolution, 

R(q) = G(-rj)b + (1 - b) , (6) 

where b is the portion from the SLR's. It is understood (or, rather, implicitly assumed) that 
the second term of light hand side, (1 — b) originates from decays of LLR's like rj, Kq , A, 
and so on 0. Table | shows the results of analysis of new data (for 8 /im target) for 
region q > 3 MeV/c using this method 0. 



III. NEW METHOD OF COULOMB WAVE FUNCTION CORRECTIONS 

We would like to propose the following new method of Coulomb wave function correction 
(with a source function p(r)) to be used instead of the Gamow factor and apply it to the 
same data as above. As usual we decompose the wave function of unlike charged bosons 
with momenta p\ and P2 into the wave function of the center-of-mass system (cm.) with 
total momentum P = |(pi + Pa) and the inner wave function with relative momentum 
q = (pi — p 2 ). This allows us to express Coulomb wave function ^(cL r ) i n terms of the 



confluent hypergeometric function F [10 



* c (q, r) = T(l - zr/)e m?/2 e iq - r/2 F(ir/; 1; iqr(l - cos 0)/2) , (7) 

where r = r\ — r 2 and the parameter r] = ma/q. Assuming factorization in the source 
functions, p(r 1 )p(r 2 ) = p R (R)p T (r) (with R = |(r x + r 2 ), / p K (R)d 3 R = 1 is assumed), we 
obtain the expression for Coulomb correction for the system of 7r + 7r~ pairs identical (modulo 
the sign) as in Refs. gflj: 
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Cc(-V) = J Pn(R)d 3 R J p r (r)d 3 r|v|/ c (q,r)| 2 (8) 

oo oo ( _ r j\n(-\m 

= G(-V) E E ^TT^ +m Hn,m)A n A* m 

n=0 m=0 11 "I" "I" 1 

= G(-i 7 )[1 + A 1o (-j7)], (9) 



where 



J(n,m) = 4vr y rfrr^^r), A n 



r(— ii] + n) 1 



r(-ZTy) (n!) 2 ' 

For the specific choice of Gaussian source distribution, p r (r) = -Cj exp(— (3 2 r 2 ), we have 



„ / i \ n+m „ 

2 / 1 \ ^ + m + 3 



/>,m) = ^M r(^^), (10) 



Using now the same method of Gaussian random numbers as in a previous Section, we can 
analyse the old and the new data on 7r + 7r~ pairs |||5j using the following formula: 

R(q)=C c (- V )b+(l-b). (12) 

Fig. ^ and Table | show our results obtained using Eq. (||) applied to old and new data with 
q > 3 MeV/c. Notice that now we can estimate (albeit with large errors) also the dimension 
of the reaction region, not available when using only Gamow factor. 



IV. STRONG INTERACTION EFFECT 



We shall consider now the possible effect of strong final state interactions, using a for- 



mulation proposed by Bowler in Ref. ||11|| . An extended formula for charged particles has 
been proposed in Ref. fll2"] . As far as our knowledge is correct, data of phase shift of p-wave 
at very small q region have not reported. Thus we consider only s-wave. The wave function 
including Coulomb effect and strong interaction effect of s-wave is described by 



*totai(q, r) = c (q, r ) + $ st(q, r) . 



(13) 



Here $ st (q, r) stands for the wave function induced by the strong interactions [I2"| and is 
given by 



4> { o'°\8)[Go(-V, qr/2)+iF (- V ; qr/2)} 



(14a) 



'(2,0) 



large r yq 



(9) exp[ i(qr/2 + r\ ln(gr) + cr ) 



(14b) 



where <t = argr(l — ir/), F and G are the regular and irregular solutions (/ = 0; s-wave) 
of the radial coordinates including the Coulomb potential ||13|| . Eq. (|14]) is the scattering 
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amplitude and the asymptotic form []. We assume that an interpolation to the internal 
region [0, r c m (1 — 2) fm] is possible as in Refs. 0,0] • (In other words, we assume a small 
contribution from the internal region in the spatial integration. In the present calculation, 
we do not introduce the cutoff factor r c .) 
The amplitude 4>q^ is decomposed as 

= \f^\0) + 1^(9), (15a) 

(£) = I e ^[exp(2^ 2) )-l] , (15b) 
%q 

/ ( i 0) W = -e J<TO [exp(2z^ 0) )-l] . (15c) 
%q 

For the phase shift of I = 2 and channel tctc scattering, we use the following parameteri- 
zation : 

1 (2) 

*° 2> w4 l + 0.5[GeVw ' ^ = - L2 ° GeV "' 

$\q) = lai 0) g, aj 0) = 1.50 GeV- 1 . (16b) 

Figs. |3| show results of our analysis of the new data. In table [I], we also show our results ob- 
tained using Eq. (|T3| ) instead of Eq. (^) in Eq. (|9|). As can be seen we observe systematically 
larger values of the interaction region with strong final state interactions included. 



V. CONCLUDING REMARKS 

We have proposed new method for the Coulomb wave function correction with mo- 
mentum resolution and applied it to the analysis of the precise data provided by Refs. 



1 Proof of Eq. ( |l3| ) with Eq. (14) is given as follows: 

e iao e iSo (cos5 ■ F + sm5 ■ G ) = e ia °[F + e iS ° sm5 ■ (G + iF )). 

Using the above equality, and summing up all the partial waves with 07 = argr(Z + 1 — irj), we have 
Eq. (|i; 

00 
1=1 

+e ir7 °e i5 ° S m5 -(G + iF )]/(qr/2) 
= * c (q,r) + $ st (q,r) . 
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Illl. Authors of Ref. || have analysed their 7r + 7r~ correlation data using Gamow factor for 
Coulomb corrections together with the random numbers method to account for final mo- 
mentum resolution. We have repeated this analysis replacing Gamow factor by the Coulomb 
wave function but following the same method for correction for the momentum resolution 
effect (cf. Eq. (|9])). As a result we were able to estimate the range of interaction for which 
(more precisely, for the root mean squared ranges of interaction r rms ) we have obtained the 
following values (for the Gaussian source function): 

_ J 9.8 ± 5.8fm for 8 fim , , , 

r rms - 2/3 ~ I 7.7 ± 4.8 fm for 1.4 /im . ^ > 

From formula including also strong final state interaction effect (cf. Eq. (|i~3"|)) we have 
obtained 



\/3 _ f 12.5 ± 4.9fm for 8 /im, 
2/3 ~ { 10.5 ± 3.9 fm for 1.4 /im. 



(18) 



Values of Eq. (^) are fairly bigger than those of Eq. (|17]). This fact means that r rms 
estimated by means of Eq. flT3"| ) depends strongly on phase shift values. The dependence is 
shown in Figs. |j. The results are also influenced by method of momentum resolution. To 
confirm Eq. (|T8|), we need to investigate different approaches in a future, for example, Ref. 
Table [Tl] shows the results concerning cutoff dependences of momentum resolution. 
The present study of 7r + 7r~ pair correlations has shown therefore that one can estimate 
the interaction region even from the 7r + 7r~ correlation data (although present data lead to 
large errors for r rms ). It can be compared with the size of the Ta nucleus, which is given by: 

(r Ta ) = 1.2 x A 1/3 

= 1.2 x (181) 1/3 = 6.8 fm . (19) 

As one can see, r rms is significantly bigger than (r Xa ). We attribute this difference to a 
physical picture shown in Fig. 0(b), i.e., to the fact that unlike-sign pions are about 60 
% emerging from the SLR's (ip, A, • • •) shown there. In a future one should apply our 
theoretical formula to other data and estimate /3's and 6's in similar reactions, and also 
consider a possibility of more direct estimation of the parameter b and its role in determining 



the source size parameter [16]. 
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FIG. 2. Results of the x 2 fits for p + Ta -> vr+vr" + X reactions with g > 3 MeV/c by Eq. (U) 
with (0): (a) 8 /im target; (6) 1.4 //m target. 
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FIG. 3. Results of the x 2 fits for p + Ta -> vr+vr" + X reactions with q > 3 MeV/c by Eq. (U) 
with (|l3|): (a) 8 /iin target; (6) 1.4 //m target. 
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FIG. 4. Phase shift dependence of 1/2/3. The horizontal axis expresses strength of a and 
a^: (a) 8 fim target; (b) 1.4 p target, TOq 2 '* and ra^ (0 < r < 1.0). 



11 



TABLES 



TABLE I. Results of the x 2 fits of R(q) for Gaussian source by Eqs. (|) and ©. 



Reaction 


Formula 


1/2/3 [fm] 


b 


X 2 /NDF 


data of Ref. § (cf. Ref. g) 


Gamow factor 


— 


— 


57.8/40 


(without momentum resolution) 


Eq. © 


2.30 ±0.88 


— 


51.0/39 


data of Ref. §§ 


Eq. @ 


— 


1.0 (fixed) 


53.6/37 


(with momentum resolution) 


Eq. (|TJ) with ® 


2.93 ± 1.03 


1.0 (fixed) 


44.7/36 




Eq. (H) with @ 


4.63 ± 0.75 


1.0 (fixed) 


44.8/36 


data of Ref. E| , 8 fim 


Eq. © 




0.43 ± 0.03 


55.0/46 




Eq. © with © 


5.63 ± 3.34 


0.53 ±0.07 


51.5/45 




Eq. (0) with ((TJ) 


7.22 ± 2.84 


0.54 ± 0.07 


51.3/45 


data of Ref. §[,1.4 fj,m 


Eq. @ 




0.51 ±0.04 


37.9/46 




Eq. (H) with © 


4.44 ± 2.80 


0.60 ± 0.07 


35.1/45 




Eq. (O) with pi) 


6.06 ± 2.23 


0.61 ±0.07 


35.1/45 



TABLE II. Momentum resolution and cutoff. 



Data 


using Monte Carlo events 


1/2/3 [fm] 


b 


X 2 /NDF 


data of Ref. || , 8 /jm 


without momentum resolution 
(central value is used) 

<Zl(T) (random) — <?L(T)| < °"l(T) 
|<Zl(T) (random) — 9l(T)| < 2<T L ( T ) 

all 


7.95 ± 3.05 

7.66 ± 2.96 
7.48 ± 2.91 
7.22 ± 2.84 


0.58 ±0.08 

0.56 ±0.08 
0.55 ±0.08 
0.54 ±0.07 


50.5/45 

50.7/45 
51.1/45 
51.3/45 


data of Ref. ||, 1.4 /um 


without momentum resolution 
(central value is used) 

l<?L(T) (random) — ^?L(T) | < C"l(T) 
|?l(t) (random) — 9l(T)| < 2<7 L(T ) 
all 


6.45 ± 2.36 

6.29 ± 2.30 
6.18 ±2.27 
6.06 ± 2.23 


0.64 ±0.08 

0.63 ± 0.08 
0.62 ± 0.08 
0.61 ±0.07 


34.7/45 

34.6/45 
35.3/45 
35.1/45 



12 



